Abstract 



We develop the theory in which rehc gravitational waves and primordial den- 
sity perturbations are generated by strong variable gravitational field of the 
early Universe. The generating mechanism is superadiabatic (parametric) am- 
plification of the zero-point quantum oscillations. The generated fields have 
specific statistical properties of squeezed vacuum quantum states. Macro- 
scopically, squeezing manifests itself in a non-stationary character of variances 
and correlation functions of the fields, periodic structures of the metric power 
spectra, and, as a consequence, in oscillatory behavior of the higher order 
multipoles C/ of the cosmic microwave background anisotropy. We start with 
the gravitational wave background and then apply the theory to primordial 
density perturbations. We derive an analytical formula for the positions of 
peaks and dips in the angular power spectrum /(/ -|- 1)C/ as a function of /. 
This formula shows that the values of I at the peak positions are ordered 
in the proportion 1:3:5:..., whereas at the dips they are ordered as 
1:2:3:.... We compare the derived positions with the actually observed 
features, and find them to be in reasonably good agreement. It appears that 
the observed structure is better described by our analytical formula based on 
the (squeezed) metric perturbations associated with the primordial density 
perturbations, rather than by the acoustic peaks reflecting the existence of 
plasma sound waves at the last scattering surface. We formulate a forecast 
for other features in the angular power spectrum, that may be detected by 
the advanced observational missions, such as MAP and PLANCK. We tenta- 
tively conclude that the observed structure is a macroscopic manifestation of 
squeezing in the primordial metric perturbations. 
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I. INTRODUCTION 



A direct search for relic gravitational waves is one of the important goals of the forthcom- 
ing gravitational wave observations (for reviews, see Refs. 0, 0). Relic gravitational 
waves are inevitably generated by strong variable gravitational field of the very early Uni- 
verse through the mechanism of superadiabatic (parametric) amplification of the zero-point 
quantum oscillations [Q. The word "superadiabatic" emphasizes the fact that this effect 
takes place over and above whatever effects occur during very slow (adiabatic) changes. 
That is, we are interested in the increase of occupation numbers, rather than in the grad- 
ual shift of energy levels. The word "parametric" emphasizes the underlying mathematical 
structure of the wave equations. It is a sufficiently quick change of a parameter of the 
oscillator, namely, variation of its properly defined frequency, that is responsible for the 
considerable increase of energy of that oscillator. 

Apparently, Schrodinger ^ was the first to notice the "alarming phenomenon" in an 
expanding universe. Specifically, Schrodinger discusses the "mutual adulteration of positive 
and negative frequency terms in the course of time". The frequency mixing means that a 
traveling wave can be amplified, with the simultaneous appearance of a "reflected" wave, 
i.e., a wave traveling in the opposite direction. After Schrodinger, out of unawareness of his 
work, this effect has been rediscovered several times. Schrodinger speaks about the mutual 
adulteration of electromagnetic waves, which would mean the generation of photons. We 
now know that the coupling of the electromagnetic field to gravity is such that the generation 
of photons is impossible, so that the alarming phenomenon does not take place. A detailed 
study of the Schrodinger paper shows that, in fact, he was operating with a variant of 
scalar electrodynamics, that is, with a scalar wave equation in an expanding universe model 
(for a discussion of this point, see Ref. P]). Then, indeed, the coupling of a scalar field 
to gravity can be chosen in such a way (minimal coupling), that the generation of scalar 
particles becomes possible. Parker 0] undertook a systematic study of the quantized version 
of the scalar wave equation in FLRW (Friedman-Lemaitre-Robertson- Walker) cosmologies. 
For a summary of the subject, see Ref. |@, 0. As for gravitational waves, there is no 
ambiguity in their coupling to gravity because the coupling follows directly from the Einstein 
equations. It was shown that the gravitational wave equation for each of the two polarization 
components is exactly the same as the equation for the minimally-coupled massless scalar 
field [0. The early studies were concerned with free test fields superimposed on a given 
space-time, whereas we are interested in fields arising in the context of perturbed Einstein 
equations (cosmological perturbations). This distinction is especially important for the issue 
of quantum normalization of the fields. 

Already at this elementary level of discussion, one can make an important observation 
that will play a crucial role in our study below. If a classical traveling wave, of any physical 
nature, is going to be strongly amplified, the resulting wave-field will form an almost standing 
wave. A traveling wave can never convert itself into a strict standing wave, because of the 
conservation of linear momentum. But the final amplitudes of the amplified left-moving and 
right-moving waves will be large and almost equal, so they interfere to form a practically 
standing wave. 

The amplification process is linear, and the final amplitude of a classical wave is propor- 
tional to the initial amplitude. If the amplitude of a classical oscillator is zero initially, the 
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oscillator will not get excited by the parametric influence. However, a quantum oscillator in 
its vacuum state does possess tiny "zero-point" quantum oscillations. One can think of these 
vacuum oscillations as the ones that are being amplified. The generation of relic gravita- 
tional waves (as well as the generation of other cosmological perturbations, discussed below) 
is a genuine quantum-gravity process, in the sense that the final result inherently contains 
all the fundamental constants h, G, and c. The gravitational energy-momentum tensor 
contains G and c, while the Planck constant h enters through the requirement of having 
initial energy ^hcu per mode of the perturbation field. The fundamental constants naturally 
combine in the Planck length Ipi = {Gh/c^Y^'^ or the Planck mass mpi = {hc/Gy^'^, but Ipi 
or nipi must stay in the numerator of the final expression, not in its denominator, so that 
the final result vanishes if h is formally sent to zero. 

Gravitational field of a FLRW universe is given by the metric 

ds'^^-c^dt'^ + a^(t)gijdx'dx^ ^a'^{r])[-dr]^ + gijdx'dx^] , (1.1) 

where the scale factor a{t) (or a{r])) is driven by matter distribution with some effective 
(in general, time-dependent) equation of state. The scale factor has the dimensionality 
of length, while r) and x^ are dimensionless. Without restricting in any way the physical 
content of the problem, one can write the perturbed gravitational field of a FLRW universe 
(for simplicity, spatially fiat) as 

ds^ = a\r))[-dr]^ + {Sij + hij)dx'dx% (1.2) 



C f°° s 1 



(1.3) 



The functions hij{r],^ have been expanded over spatial Fourier harmonics e^'"'', where 
n is a constant (time-independent) wave vector. The wave number, n, is related to n by 

n = {SijTtn^Y^'^ . The wave number n defines the wavelength measured in units of laboratory 
standards (so to say, in centimeters) by A = l-najn. Using the Fourier expansion, we are 

able to reduce the perturbed dynamical problem to the evolution of mode functions hnij]) 
for each mode n. Two polarization tensors Pjj(n), s = 1,2 have different forms depending 
on whether they represent gravitational waves, rotational perturbations, or density pertur- 
bations. If (n/n, 1, m) are three unit and mutually orthogonal (spatial) vectors, then we 
have for gravitational waves, 

1 2 S ; ■ S ■ 

Pij = ~ mjini, Pij = hrrij + Ijirii, PijS'^ = 0, Pijn^ = 0, 
for rotational perturbations, 

1 1 2 1 s ,. s ■ ■ 

Pij = -{krij + IjUi), Pij = -{mirij + mjUi), Pijd'^ = 0, PijTi'n^ = 0, 
and for density perturbations, 

Pij = y gdij, Pij = + 
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In all three cases, Pjj(n) obey 

Pij{n)p '^{n) = 2Sss', Pij{-n) = P^j(n). 

In general relativity, rotational and density perturbations can only exist if they are 
supported by the corresponding perturbations of matter. Their propagation speeds depend 
on the properties of matter and can range from zero to the speed of light, c. For instance, the 
propagation speed of density perturbations in the radiation-dominated fluid is c/ -\/3; and 
it was very close to c if the very early universe was driven by a scalar field However, 
in alternative theories of gravity, solutions with the polarization structure of rotational and 
density perturbations can exist even in the absence of matter fields, in which case the metric 
perturbations represent gravitational waves with new polarization states, in addition to 
the usual gravity- wave polarization states of general relativity If one concentrates on 
metric perturbations alone, temporarily leaving aside the accompanying perturbations of 
matter variables, then all three types of cosmological perturbations in general relativity can 
be thought of as gravitational waves, even though some of them have unusual polarization 
states and unusual propagation speeds. There is no wonder that the dynamical equations for 
cosmological rotational and density perturbations are similar to, and sometimes exactly the 
same as, equations for cosmological gravitational waves. The common "master equation", 
whose solutions allow one to derive all the metric components along with all the matter 



perturbations (when they are present), has the universal form |jTT 

/" + / 



n4 - 



:i.4) 



where / := d/dr], q is a function of rj and is interpreted as the propagation speed of the 
perturbation, and W{7]) is a function of a{r]) and its derivatives. For density perturbations 
in a perfect fluid with the fixed equation of state p = we, ci is a constant. In the case of 
gravitational waves, cf /c^ = 1 and W{ri) = a" /a M. One can view Eq. (|1.4|) as the equation 



of an oscillator with variable frequency (the term in square brackets), or as the Schrodinger 
equation of a particle moving in the presence of a potential barrier W{7]) (while remembering 
that ?7 is a time coordinate rather than a spatial coordinate). In what follows, we will be 
discussing gravitational waves and density perturbations. 

For a classical gravitational field, the quantities c^, in Eq. ( |1.3|) are arbitrary complex- 
conjugate numbers. The constant C can be incorporated into them. In the quantized version, 

the quantities c^, are annihilation and creation operators satisfying the conditions 

[c„, cj = 6s'sS%n - m) , c„|0) = , (1.5) 

where |0) (for each n and s) is the fixed initial vacuum state defined at some r/o in the very 
distant past, long before the superadiabatic regime for the given mode has started. In that 

s 

early era, the mode functions hniv) behaved as oc e"™'', so that each mode n represented 
a strict traveling wave propagating in the direction of n. The normalization constant C is 
VTOvr/pi for gravitational waves, and \/2A^lpi for density perturbations. 
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A detailed study shows [0 that the quantum-mechanical Schrodinger evolution brings 
the initial vacuum state of cosmological perturbations into the final multi-quantum state 
known as squeezed vacuum state. It is the variance of phase that is being strongly di- 
minished (squeezed), while the mean number of quanta and its variance are being strongly 
increased. A squeezed vacuum state is conveniently characterized by the squeeze parameter 
r. The squeeze parameter grows from r = in the vacuum state up to r ^ 1 by the end of 
the amplifying superadiabatic regime. The mean number of quanta in a 2-mode squeezed 
vacuum state is (A^) = 2sinh^r. Squeezed vacuum states possess specific statistical prop- 
erties. In particular, the generated field, viewed as a random field, obeys the statistics of a 
Gaussian non-stationary process. The non-stationarity means that the variance of the field 
is an explicit oscillatory function of time, and the two-time correlation function depends on 
individual moments of time, not only on the time difference. The calculation of quantum- 
mechanical expectation values and correlation functions provides the link between quantum 
mechanics and macroscopic physics. 

Using the representation (|1.3|) and definitions above, one finds the variance of metric 
perturbations: 

^oo ^ s dr; 

m,{v,^)hHv,^m = ^J E \hn{v)\'-- (1.6) 

/TT JO s=l,2 ^ 

The quantity 

h\n,v) = ^y E \Uv)\' (1.7) 

s=l,2 

gives the mean-square value of the gravitational field perturbations in a logarithmic interval 
of n and is called the (dimensionless) power spectrum. In the case of gravitational waves, it 
is relatively easy to evolve the mode functions up to the present era, and to find that 

/i^ (n, r]) oc sin^ [n(r7 — ?7e)], (1.8) 

where rje is a constant discussed below. The explicit time-dependence of the power spectrum 
is a consequence of squeezing and can be also viewed as a refiection of the standing-wave 
pattern of the generated field. For every fixed moment of time (for instance, today) the 
power spectrum contains many maxima and zeros at certain wave numbers, even though 

s 

the spectrum was perfectly smooth before amplification, i.e., when the mode functions hniv) 
behaved as oc e"*"''. As soon as the amplifying process takes place, the increase of the 
mean number of quanta, squeezing, non-stationarity, formation of standing wave pattern 
and oscillatory features in the power spectrum, - are all the different facets of the same 
phenomenon. 

The relative spacing of zeros is very dense at laboratory scales (large n's), but is quite 
sparse at cosmological scales (small n's). Specifically, the spectrum contains about 10^° 
zeros in the interval from 100 Hz to 200 Hz, but only a dozen of zeros in the interval from 
1000 Mpc to 2000 Mpc. The oscillatory time-dependence ( |1.8| ) is known in advance, and this 
information would certainly help, in a very narrow-band gravitational wave detector, to find 
the signal against the instrumental noise, and to provide evidence for the primordial origin 
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of the detected gravitational wave background [Q]. However, in a broad-band detector, there 
are too many zeros together, and the non-stationary process is practically indistinguishable 
from the stationary process of the same power density. In a recent paper, Allen, Flanagan, 
and Papa [jl3[ agree that the non-stationarity in the relic background is present, but they 
argue that this feature can hardly give any advantage in practice. Without presently being 
able to offer a realistic scheme of exploiting the non-stationarity at small scales, we shift 
our attention to cosmological scales, where the spacing of zeros is sparse. The natural place 
to look for the consequences of squeezing is the distribution of multipoles of the cosmic 
microwave background (CMB) anisotropics M . 

Relic gravitational waves should be important at lower multipoles I of the CMB 
anisotropics, but they are not expected to give significant contribution at £ ~ 200 and 
higher. It is the density perturbations that are expected to be primarily responsible for 
the behavior of the C^'s in the latter range. However, we study in great detail the sim- 
pler case of gravitational waves in order to resolve a number of principal issues and to get 
a guidance for the analysis of the technically more difficult case of density perturbations. 
The phenomenon of squeezing is universal, and if the primordial density perturbations have 
a quantum- mechanical origin (which, we believe, is likely to be true), then many features 
must be common with the case of relic gravitational waves. Some differences arise at the 
late stages of evolution (in particular, they explain why the gravitational wave [g.w.) contri- 
bution is subdominant at ^ ~ 200 and higher), but we take them into account. Developing 
the conjecture of Ref. |T^, we argue that it is the modulated structure of the metric power 



spectrum (which is caused, on fundamental level, by squeezing) that is responsible for the 
downturn of the rising function £{£+ l)Ce at the peak, and for the appearance of subsequent 
peaks and dips, a few of which have been recently observed [IG]. On the ground of our simple 
analytical treatment we make a forecast for the positions of further peaks and dips that may 
be observed by future missions, such as MAP and PLANCK. In general, our forecast agrees 
with that of Ref. |T^, made on the grounds of numerical codes, but discrepancies become 
significant somewhere around the 4-th expected peak. 

The structure of the paper and its conclusions are as follows. We start, in Sec. ||, with the 
gravitational wave solutions in the present universe, that is, at the matter-dominated stage. 

s 

We consider the general solution for the time dependent mode functions hn{v)y regardless of 
whether a given particular solution is likely to emerge from the very early universe or not. 
In general, a given mode n is neither a traveling wave nor a standing wave. We formulate 
conditions under which a given mode is a strict traveling wave or a strict standing wave. 
These conditions are constraints on the (Fourier) coefficients An, Bn in front of two linearly 
independent solutions for the mode functions. Then, we explore the issue of stationary ver- 
sus non-stationary variance. We demonstrate that oscillations in the power spectrum are 
most pronounced when the modes are standing waves, and they disappear when the modes 
are traveling waves. Thus, the often made (incorrect) statement that the relic gravitational 
wave background should be stationary is equivalent to the assumption that it is being formed 
by traveling waves. We return to this issue later on (in Sec. fV\j and show that, whether the 
non-stationarity on small scales is measurable or not, the very assumption of the stationary 
gravitational wave background is in conflict with some other cosmological considerations. In 
Sec. |T|, we present a "physical" model for the gravitational wave background. We call it 
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"physical" because we evolve the field through all the three relevant stages of cosmological 
evolution - initial, radiation-dominated, and matter-dominated. In this way we distinguish 
the "physical" model from the "alternative" model, which postulates that the gravitational 
wave background is stationary, irrespective of its physical origin. The waves generated in 
the physical model are squeezed (standing). We demonstrate that the evolution of stand- 
ing waves through the effective barrier a" /a at the matter-dominated stage results in the 
appearance of an oscillatory behavior of the Fourier coefficients An, Bn as functions of the 
wave number n. We later show (in Sec. 0) that this oscillatory behavior of An, -B„ is the 
origin of oscillations in the multipole moment distribution Ce as a function of i. 

Section ^ compares the physical and alternative gravitational wave backgrounds. We 
introduce the notion of a fair comparison, which requires that today's band-powers of the 
two backgrounds be equal at all scales. By evolving the corresponding solutions backwards 
in time, we show that the alternative background would have had too much power in long 
waves at the era of last scattering of the CMB radiation. In terms of "growing" and "de- 
caying" solutions, this means that the amplitude of the decaying solution in the alternative 
background becomes dangerously large when one returns deeper and deeper in the past. 
The further evolution backwards in time would have destroyed our sacred belief that the 
Universe was homogeneous and isotropic (up to small perturbations) at the time of the 
primordial nucleosynthesis and its past. Most importantly for our study, we demonstrate 
that the alternative background does not produce oscillations in the Cg multipoles. This 
shows that the squeezing is observationally distinguishable, even if, at the present level of 
observational capabilities, it is better to search on large scales rather than on small scales. 

Section [V| presents the results of numerical calculations of the Cgs caused by relic grav- 
itational waves. We show that our analytical formula for the positions of peaks and dips is 
in a fairly good agreement with numerical calculations. This analysis demonstrates that, at 
least in the case of gravitational waves, the Ce oscillations are produced by modulations in 
the power spectrum of metric perturbations, and not by acoustic waves at the last scattering 
surface, simply because there are no matter perturbations at all. 

In Section we turn to the primordial density perturbations. The evolution of density 
perturbations through the initial and radiation-dominated stages is almost identical to the 
evolution of gravitational waves. We show that the Fourier coefficients An, Bn of metric per- 
turbations associated with density perturbations develop a periodic structure, as functions 
of n, in the course of transition from the radiation-dominated phase to the matter-dominated 
phase. In a manner similar to the gravitational wave case, we derive the expected positions 
of peaks and dips in the Q distribution. Because of the damping, features beyond the sec- 
ond peak may not be easily discernible |]TH|, |[TD|, [^, but they seem to be less likely to be 
washed out if they are produced by modulations in the metric power spectrum rather than by 
modulations in the plasma matter power spectrum. We show that the peak positions should 
obey the rule 1:3:5:7..., whereas the dip positions should be ordered as 1:2:3:4.... We demon- 
strate that the observed positions agree better with our analytical formula than with the 
concept of "acoustic peaks" . The detection of late features may be especially interesting as 
our analytical forecast is somewhat different from what follows from conventional numerical 
codes. We tentatively conclude that the observed structures in the angular power spectrum 
i{£ + l)Ce are macroscopic manifestations of squeezing in gravitational field perturbations. 
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II. PROPERTIES OF GRAVITATIONAL WAVE SOLUTIONS 



The perturbed Einstein equations for gravitational waves can be reduced to the "master 
equation" 



s'l s 



n 







(2.1) 



where the functions /^nl^?) are related to the mode functions hn{v) by 

s s 



(2.2) 



We suppress the polarization index s when it causes no ambiguity. 

The scale factor a{7]) at the matter-dominated stage, governed by whatever matter with 
the effective equation of state p = 0, behaves as a{r]) o:rf. It is convenient to write 0(77) in 
the explicit form 



a{r]) = 2lH{r] - Vrnf 



(2.3) 



where In is the Hubble radius today {In = c/Hq, where Hq is the present value of the 
Hubble parameter) and rjm is a constant explained below. The moment of time "today" 
(in cosmological sense) is labeled hj rj = rjR (the subscript R denoting "reception"). It is 
convenient to choose 



VR-Vn 



1. 



(2.4) 



With this convention, a^rjn) = 21h, and the wave, of any physical nature, whose wavelength A 
today is equal to today's Hubble radius, carries the constant wavenumber uh = 47r. Longer 
waves have smaller n's and shorter waves have larger n's, according to the relationship 
n = AttIh/X. For example, the ground-based gravitational wave detectors are most sensitive 
to frequencies around 30 Hz - 3000 Hz. The corresponding wavelengths have wavenumbers 
n somewhere in the interval 10^'^ — 10^^. 

For the scale factor (p.3|), Eq. ( p.l| ) is easily solved to yield 



2 



{A^ - Bn) Hi%{y) + (A„ + Bn) H'^/M 



(2), 



where 



y = n{r]- r]m) 



(2.5) 



(2.6) 



and J±3/2{y) and Hy2\y) are Bessel and Hankel functions, respectively ||2T|. We will also 

be using spherical Bessel functions ji{y) = ^7r/2?/Jj+i/2(l/)- The (Fourier) coefficients An 
and Bn are, so far, arbitrary complex numbers. In the limit y ^ 00, 
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e and ji{y) 



-cosy, ]-2{y) ~ --siny. 

y y 

(2.7) 



Thus, at late times, Uniji) is a combination of sine and cosine functions of time: 

1 



{An COS y - iBn sin y) 



(2.8) 



A. Standing versus traveling waves 



We now consider a classical field ( |1.3|) . The n-th mode of the field is given by the real 
function 



where a complex number Cn can be conveniently presented in its polar form: 



(2.9) 



(2.10) 



The gravitational wave solutions at late times are given by Eq. ( |2.8| ). Thus, we have 



a{r]) 



^ [ {An cos y - iBn sin y) e^("-^+'^-n) 
+ {A*n cos y + iBl sin y) e-^("-^+<^-" ) 



(2.11) 



If the coefficients An-, Bn are arbitrary, Eq. (|2.11| ) is neither a traveling wave nor a standing 
wave. A traveling wave is characterized by two real numbers, namely, an amplitude A and 
a phase 0; its general form is y4sin(±n?7 + n ■ x + 0). The minus/plus sign describes a 
wave traveling in the positive/negative direction defined by the fixed vector n. A standing 
wave is characterized by three real numbers, viz., an amplitude A and two phases, and its 
general form is Asai{nri + 0i) sin(n • x + 02)- Different choices of these free parameters are 
responsible for concrete space-time patterns, but they do not change the wave classification. 
A little investigation shows that Eq. ( p.ll| ) describes a traveling wave if and only if 



Ar 



±Bn. 



This constraint can also be written as 



Ar, 



pne 



i(j>„ 



and Br, 



i4>n 



(2.12) 



(2.13) 



where pn and 0„ are two arbitrary real numbers. We will refer to this constraint as the 
traveling wave condition. On the other hand, Eq. ( |2.11| ) describes a standing wave if and 
only if 
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An = pA„e''^" and 5„ = ±ipB^e'' 



(2.14) 



where pa„, Pb„, 0n are three arbitrary real numbers. We will refer to this constraint as 
the standing wave condition. Two special cases of the standing wave condition are when 
either y4„ = or i?„ = 0. This classification of traveling and standing waves was based on 
the regime when a given mode satisfies the requirement y ^ 1 (short-wave regime), but it 
can now be applied at earlier times too, when this requirement is not satisfied (long-wave 
regime) . 



B. Stationary versus non-stationary variance 



The variance of a quantized field ( |1.3| ) is defined by Eq. ( |1.6| ). The essential part of 
this expression is the power spectrum given by Eq. ( |1.7|) . At the matter-dominated stage, 
the general solution for the gravitational wave mode functions is represented by Eq. ( p. 51) . 
One can now find the power spectrum. We use spherical Bessel functions and replace the 
summation over s by the multiplication factor 2. Then, the general expression for the power 
spectrum is 



Anrjiiy) + \Br 



\^£,iy) + 2lm{A:B„ 



^3i{y)j-2{y) 



(2.15) 



Using the constraints ( p.l3| ) and (|2.14|) , one can now specialize the power spectrum to 
traveling and standing wave cases. In the traveling wave case one obtains 



and in the standing wave case. 



2C'nYp\ 



Jiiy)+£2iy) 



(2.16) 



2C^nV 



[pAjiiy) ± PbJ~2 



(2.17) 



When considering the waves that are shorter than the Hubble radius, ^ 1, one can use 
asymptotic formulas (|2.7|) for the Bessel functions. One can also replace a{ri) with 21 h, 
since the scale factor does not practically change during any reasonable observation time. 
Then, in the traveling wave case, the oscillations of the power spectrum fully disappear, 
as the oscillating terms cos^ y and sin^ y combine to 1. However, the oscillations are most 
pronounced in the standing wave case, as Eq. ( p.l7|) exhibits a periodic structure 



27C% 



[pA„ COS y ±Pb„ sin yf , ?/ > 1. 



(2.18) 



In other words, the power spectrum of traveling waves is stationary, whereas the power spec- 
trum of standing waves is non-stationary. These two classes of power spectra are subcases 
of the general situation in which An = PAn^^'^^"^ = P_B„e*"^^", and correspondingly. 
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Pa„ cos^ y + sin^ y + 2pa„Pb„ sin(0B,^ - sinycosy , y > 1 



(2.19) 



At a fixed moment of time, for instance today, the power spectrum (|2.18| ) reduces to 



[pA„cosn ± pB„sinnf , n > 1. 



(2.20) 



As one can see, tlie power spectrum of standing waves contains, quite generically, many 
maxima and zeros at certain n's due to the oscillatory factors cosn, sinn. The coefficients 
PAn, PBn are still arbitrary and could, in principle, be smooth functions of n. However, 
as we will show in the next section, the preceding evolution of standing waves through 
the transition from the radiation-dominated era to the matter- dominated era, gives rise 
to additional oscillations, namely to oscillations in the coefficients pA„, Pb„ themselves, as 
functions of n. 



III. THE PHYSICAL MODEL FOR RELIC GRAVITATIONAL WAVE 

BACKGROUND 

A. The behavior of the scale factor: Pump field 

The matter-dominated (m) era with the scale factor ( |2.^ ) was preceded by the radiation- 
dominated (e) era with the scale factor a{r]) oc t]. To simplify the analysis, and with- 
out any essential loss of generality, we assume that the transition from e era to m era 
was instantaneous and took place at some t] = 772. The redshift of the transition is Zeq-. 
(^{Vr)/(^{V2) = 1 + Zeq. It is bclieved that Zeq is somewhere near 6 x 10^. In its turn, the 
radiation-dominated era was preceded by the initial (i) era of expansion, whose nature and 
scale factor are, strictly speaking, unknown. To simplify the analysis, and since the wave 
equations admit simple exact solutions in case of power-law scale factors, we assume that 
the i era, similar to the e and m eras, was also described by a power-law scale factor. We 
parametrize the i era by a{ri) oc Ir^l^"^^ (compare with [Q). The transition from i era to e 
era takes place at some rj = rji and at redshift zf. aijjR) / airji) = 1 + Zi. Further analysis 
shows (see below) that in order to get the right amplitude of the generated perturbations, 
the numerical value of Zi should be somewhere near lO^*'. 

We now write the full evolution of the growing scale factor explicitly: 

a{r]) = Qv\''~'', V<Vi, Vi<0, /3 < -1 , (3.1) 
a{r]) = loae{ri - r]e), r]i < V < V2 , (3.2) 

a{r]) = 2lH{r] - ?7„)^, t]2 < i]. (3.3) 

The continuous joining of 0(77) and a'{ri) at the transition points fully determines all the 
participating constants in terms of In, Zi, z^q and (3. Concretely, 
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rjR = l + r]r, 



Vn 



eq 



1 + Zeq ' 



V2 



1+/3 



eq 



1 + Zi 



eq 



1 yrr^ 



2^ 1 + Zi ' 



and 



lo^^e 



41 



H 



- 1 



I 



1 + (1+^.)-^ ■ 



a/1 + ^eq 

The case /? = — 2 is known as the de Sitter inflation. In this particular case, 

a/1 + ^eq 



(3.4) 



(3.5) 



lo 



H- 



1 + • 



(3.6) 



For the CMB calculations we will also need the redshift Zdec of the last scattering surface 
f] = Ve (with the subscript E denoting "emission"), where the CMB photons have decoupled 
from rest of the matter: a{r]^ / a{r]E) = 1 + z^ec- The numerical value of z^ec is somewhere 
near 1000. The time of decoupling vjE is 



1 



1 



Ve 



Zde 



2^/1+^ 



eq 



a/1 + Zeq 
2il + Zi) 



(3.7) 



All the formulas above are exact and we will be using them often, but surely one can get 
an excellent approximation by neglecting 1 in comparison with ^j, Zf.q and Zdec- 



B. Squeezed gravitational waves 

As soon as the scale factor and initial conditions for the mode functions are strictly 
defined, the coefficients An,Bn in Eq. ( p.5|) are strictly calculable p2| , p^ . The general 
solution for yU„ at the e stage is 

finiv) = ^le-'^^"-"^) + Eae'"^^-''^), (3.8) 

but the preceding evolution allows one to specify Bi and i?2- The waves subject to ampli- 
fication at the i stage have the wave numbers n satisfying the condition ?T,|r7i| ^ 27r|l + 
that is, 

n 1 + Zj 



eq 



For these waves, one finds 



n^/T+Zeq 



B,^-B,^F{/3)1 \^^'' ^B, (3.9) 



where 
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_gi(xo+7r/3/2)_ 



22{/3+l)r(/? + 3/2) COS pTl 

and xq = nrjo p2|. Note that in the particular case P = —2, 



\B\ 



4(1 + Zi)^ 



-2. 



(3.10) 



(3.11) 



The fact that Bi ^ —B2 demonstrates that the gravitational wave modes n are (almost) 
standing waves at the e stage. These standing waves are encountering the nonzero barrier 
a" /a at the m stage. To find the coefficients An, Bn in Eq. (p.5|) one needs to join contin- 
uously the finiv) and fi'niv) at the transition point t] = rj2- This calculation allows one to 
find the coefficients An, Bn, but it also shows (as expected) that the standing-wave character 
of the field at the e stage leads to the appearance of oscillations in the coefficients An, Bn 
of the field at the m stage. This is a general phenomenon which we will also discuss in 
connection with density perturbations. Explicitly, 



A — ■ 1^ 



B 



Br, 



2%i 

Tf B 

2%I 



(81/2 - 1) sin y2 + 4?/2 cos y2 + sin 3?/2 
[%yl - 1) cos 7/2 - 4?/2 sin ?/2 + cos 3?/2 



(3.12a) 
(3.12b) 



where 



y2 = n{r]2 - r]e 



n 



and rir, 



2jl + Zeq. 



The numerical value of is about 160. This corresponds to wavelengths that are 15 times 
shorter than l^- Clearly, the coefficients ( 3.12|) satisfy the standing wave condition {^J^. It 
is easy to check that had one artificially chosen traveling waves at the e stage, by assuming 
that either Si = or _B2 = 0, the oscillations in An, Bn would have been suppressed. 

The asymptotic expressions for An, Bn are as follows. For relatively short waves, i.e., 
?/2 ^ 1, n ^ rim, one has 



An^ -i2^7v/2Bsmy2, 5„ ^ -2^7r/25 cos?/2. (3.13) 
Using these coefficients in Eq. ( |2.8|) , one finds the finiv) at the m stage, for n 3> rim- 

IJ'niv) ~ -i2B sin[n{r] - T]e)]. 

Of course, this function is simply a continuation of the relevant standing wave solution 
{Bi f« —B2 ~ B) at the e stage, Eq. (|3.8D, to the m stage. Indeed, the height of the barrier 
at the m stage is 



a 
a 



V=V2 



1 2 

2 '^m) 



SO the waves with n ^ Um stay above the barrier and experience no changes. However, the 
waves with n ^ rim and n <^ are affected by the barrier. For relatively long waves, i.e., 
n <^ rim, one has 
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An ~ 



2 ) 



(3.14) 



2v^ ' " 45v^ 
such that \Bn\ ^ |^n|- The formulas above are in full agreement with if one takes 



into account the change of notations: An = y7r/2Cin, Bn = i^Ti/2C2n- In particular, the 
long- wavelength part of the power spectrum at the time of decoupling rj = tje is given by 



4/j 



h'^{n, Ve) ^ --^{1 + Zeq)n^\B\^, n < ridec = Vl + Zde 

TT If 



(3.15) 



This part of the spectrum is primordial, in the sense that it has not changed since the 
beginning of amplification. In the particular case (3 = —2, and using Eq. ( p.ll| ), we obtain 
the flat (independent of n) primordial spectrum 



TT 11(1 + Zeg) ' 



/3 



n < Ude 



Zde 



(3.16) 



In preparation for the discussion of CMBR anisotropies, we show in Fig. |l] the numeri- 
cally calculated spectrum h?{n,r]E), including the beginning of its oscillations. We use the 
notations 



x = y2 



n 



b = 



eq 



nde 



VI + Zde 



The substitution of An, Bn given by Eq. ( |3.12| ) into the general expression ( [^.15| ) produces 
the exact power spectrum 



/2 



^pi 



'l + ZdecWWg\x,h) 



(3.17) 



H 



where 



g'^{x,h) = [pi{x)ji{bx) - p2{x)j^2ibx)] , x= — , b 



n 



Ude 



PlW = — 



[8x'^ — 1) sin X + 4x cos x + sin 3x 



P2[X) = — 
X^ 



(8x^ — 1) cosx — 4xsinx + cos3x 



Note that the rigorously evolved mode functions single out the lower sign in the general 
formula ( p.l7| ) for standing waves. Fig. |l| shows the function g'^{x, b) for 6 = 5. 

The positions of maxima and zeros of the power spectrum are well approximated by the 
short- wavelength limit of g'^{x, b): 



9\x,b) 



64 



sin^[(6 — l)x]. 
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FIG. 1. Plot of b = 5) versus x. 
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In the case of 6 = 5, the maxima and zeros of the function sin^ 4x are ordered as follows: 



Maxima: ^r'' = ^{2k + l), x] 



max 
k 



a;r"(2A;+l), A; = 0, 1, 2, 3, 



(3.18) 



Zeros: xr" 



-2fc, 



mm 



(3.19) 



The analytical forecast (|3.18| ), (|3.19| ) shifts the first few features to the left (smaller x) 
as compared with the more accurate numerical calculation, but the forecast becomes pro- 
gressively more accurate for later features. For example, the analytical formula predicts 



XV 



0.79, 



XV 



1.18 



) ^8 ~ 
min _ Q_Qi 



places these features at 



6.28, 



Xo 



6.68, whereas the numerical calculation 



, X^ 



1.21, 



Xq 



6.29, 



Xo 



6.68. The zero-th 



maximum x, 







would be placed, according to the analytical formula (|3.18| ), at x^""^ = 0.393. 



The crucial rule, however, is that the positions of minima, starting from x^***^, are ordered 
in the proportion 1:2:3:4.., whereas the positions of maxima, starting from the zero-th 
"would-be" maximum, are ordered in proportion 1:3:5:7.... We will see in Sec. |V| how 
these features are being reflected in the oscillations of the Ci multipoles. 



IV. DISTINGUISHING BETWEEN THE PRESENCE AND THE ABSENCE OF 

SQUEEZING 

As demonstrated above, the quantum-mechanically generated (squeezed) gravitational 
waves form a non-stationary background, whose power spectrum is fully determined by fun- 
damental constants and parameters of the gravitational pump field. In this section we show 
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that the hypothesis whereby the gravitational wave background is postulated to be station- 
ary, and therefore non-squeezed, is in conflict with some other cosmological data. Since we 
concentrate on stationarity versus non-stationarity, the comparison of the two backgrounds 
should be fair, in the sense that their today's (for r] = rj^) band-powers should be assumed 
equal. However, as we will show below, when one returns back in time to, say, the decoupling 
era, the alternative background is bound to have too much power in long- wavelength per- 
turbations. The further extrapolation back in time destroys the usual (and partially tested) 
assumption that the cosmological perturbations remain small all the way down to the nu- 
cleosynthesis era and beyond. In Sec. |V], we will also show that the alternative background 
does not produce oscillations in the Ci multipoles, whereas the physical background does. 

The power spectrum of the physical background is given by the general expression (p.l5|) 
where An,Bn are determined by Eq. ( p.l2| ). The substitution of Eq. ( |3.12| ) in Eq. ( p.l5| ) 
results in formula (|2.17|) (with the lower sign) where pa„ , Pb„ are absolute values of the 
coefficients ( |3.12| ). On the other hand, the general form of the power spectrum for the 
alternative (stationary) background is given by Eq. ( p.l6| ), where pn should be found from 
the fair comparison: 



jiH+£2( 



n] 



[PAnJiin) - PB„J-2{n)Y 



(4.1) 



We will find p„ using asymptotic expressions ( p.l3| ), ( p.l4| ) in different intervals of n. In the 
region n S> rim, formula (|3.13| ) yields 



PA„ 



'2ti\B\ 



n 



sm 



PBr, 



2tx\B\ 



n 



cos ■ 



(4.2) 



In the region n <^nm formula ( |3.14|) yields to 

3v/7r rim 



71/ rn. 



(4.3) 



Let us start from n 3> rim- Using asymptotic expressions ( |2.7|) for the spherical Bessel 
functions, and replacing sin^n with 1/2 in the right-hand-side of ( |4.1| ), one finds 



(4.4) 



Now turn to the interval rim ^ n ^ 1. The p^^, pb„ are now given by (O). One can 
still use the asymptotic formulas ( |2.7|) , but the second term in the right-hand-side of ( [4.1|) 



is much smaller than the first one and therefore can be neglected. Then, one derives 



pl^^lBl""'-^, n^>n>l. 
Id ""^ 



(4.5) 



Finally, in the region n <C 1, one uses asymptotic formulas ji{n) n/3, j-2{n) ~ l/n"^. 
In either side of ( [4.1|) , the first term is smaller than the second term, and can be neglected. 
Then, one obtains 



71 



pl^-\B\^nln\ n<l 



(4.6) 
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These formulas give a piece-wise representation for the smooth alternative spectrum which 
today has, in all intervals of n, approximately the same power as the physical spectrum does. 

Since the coefficients p^,^, p„ are fully determined, one can now derive the forms of 
the two power spectra at other times. We will compare the two backgrounds at the time of 
decoupling rj = rjE- One needs to consider formulas (|2.16|) and ( |2.17| ), where a{ri) = a{rjE) 
and y = i/e, 



Ve = n{r]E - r]„ 



n 



ride 



n 



(4.7) 



Clearly, in the band of sufficiently short waves, n ^ a/1 + Zdec-, both spectra increase power 
in the same proportion, simply as a result of changing of 0(77) from a{rifj) = 21h to the 
smaller value aijjE) = 21^/(1 + Zdec)- So, in this range of n the ratio of powers in the two 
backgrounds is 1: 



h'^{n,r]E)\ait 
h'^{n, r]E)\phys 



1, n > VI + Zde 



(4.8) 



However, this ratio is significantly larger than 1 in longer waves. This is seen from the 
general formula 



pI 



h^{n,r]E)\ait 

h'^{n, r]E)\phys [pAji{n/ndec) - pBj~2{n/ndec)\ 



(4.9) 



applied to longer waves. Indeed, in the interval 1 <^ n <^ a/1 + ^dec (and, hence, n -C rim) 
one uses the small argument approximation for the Bessel functions, neglects second terms 
in numerator and denominator, and takes from Eq. ([4.5|) . This calculation results in 



h'^{n,r]E)\ait 
h'^{n,riE)\phys 



9 
2' 



1 + Zde 



K < VI + Zde 



(4.10) 



This ratio is comparable with 1 only at n ~ VI + Zdec where ([4.10|) goes over into (^^ ). But 
the ratio ( ^.10| ) is much larger than 1 for smaller n's. As for the region n -C 1, one applies 
the same approximations as in the previous case, but takes from Eq. ( [4.6|) . This gives 
the result 



h\n, riE)\ 



alt 



h'^{n,T]E)\ 



phys 



{1 + Zdee?, n-^l. 



(4.11) 



This ratio is universally (independently of n in this region) much larger than 1. 

It is easy to understand these results. The parts of the general solution ( |2.5| ) with ji (y) 
and j-2{y) are called, respectively, the growing and decaying solutions. This classification 
refiects their different behavior in the small argument approximation for the Bessel functions. 
The growing and decaying solutions are necessarily present in the power spectra of both 
backgrounds. However, the stationary background contains jf{y) and ^^2(1/) always with 
equal coefficients, whereas the non-stationary background contains ji{y) and j-2{y) with 



equal (up to oscillations) coefficients (|4.2|) only for sufficiently short waves. In longer waves, 
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the decaying solution of the stationary background becomes progressively more important 
when one goes back in time. Since the physical background is expected to have enough power 
to produce the lower order anisotropy in the CMB temperature at the level 6T /T ~ 10~^, 
the alternative background would have produced this anisotropy at the unacceptably high 
level 5T/T ~ 3 x 10^^ 



V. MICROWAVE BACKGROUND ANISOTROPIES CAUSED BY RELIC 

GRAVITATIONAL WAVES 

The key element in formula for the temperature variation 5T /T seen in a given direction 



25| is the ?7-time derivative of metric perturbations evaluated along the CMB photon's 



path between the event of reception (R) and the event of emission (E): 

5T , . 



T 



The upper limit of integration is 



1 
2 







Jo 


dr] 



dw. 



(5.1) 



- path 



Wl='qR- 7]E 



1 - 



and the integration is performed along the path r] = rj^ — w, x = ew. In the case of 
density perturbations, the integral in formula ( |5.1| ) should be augmented by the additive 
term representing initial conditions: an intrinsic variation of temperature at rj = rjE and 
a possible velocity of the last scattering electrons with respect to the chosen coordinate 
system, which is synchronous and comoving with the perturbed, gravitationally dominant 
pressureless matter, possibly - cold dark matter (CDM). 

Similar to the perturbation field hij itself, the temperature variation 6T/T is also a 
quantum-mechanical operator. To establish contact with macroscopic physics, we need to 
calculate the correlation function 







5T 5T 
— (ei) — (ea) 







We use the mode functions (|2.5|) and the normalization constant C = ^/TQttIpi. Then, it can 
be shown [^] that the correlation function takes the elegant form 







6T 6T , 
— (ei) — (62, 



0) =/|,X]ir,P/(cos5), 



(5.2) 



1=2 



where P;(cos5) are Legendre polynomials for the separation angle 6 between the unit vectors 
ei and 62, and 

Ki = (2/ + l)(/-l)/(/ + l)(/ + 2) Fi, 



where 
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Fi 



n 



Ji+i/2{nw) 



[nw 



)5/2 



fniVR - W)dw 



and 



1 (f^nV 

\/2n\ a J 



v=vr-w 

We will also be using the multipole moments Ci defined by 

21 + 1 



(5.3) 



(5.4) 



4:71 



-G 



The central quantity for the calculation of Fi is the function 



a 



n 



--V^(^n ^5/2(2/) +iBnJ-5/2iy))- 



(5.5) 



To get more insight into Fi we introduce the two new functions: 'ip±i, defined by the respective 
integrals: 



"'"I ^^Ji+i/2{x) J±5/2{n - x) 



^5/2 (j^ _ 3;^3/2 

Then, the general expression for Fi takes the form 



Fi 



Jo 



dn n 



AnVri + \BnVr-i + 2im(A:s„)^i^_ 



(5.6) 



(5.7) 



This expression simplifies when the coefficients A^, represent traveling (|2.12|) or standing 
(2.14) waves. For a stationary background one gets 



Ul Jo 



while for the physical nonstationary background one gets 



Jo 



(5.8) 



(5.9) 



These formulas explain the different behavior of the multipole moments Ci in the two cases. 
We demonstrate this with the help of numerical calculations. 

In Fig. ^ we show by a solid line the graph of the function l[l + l)Ci calculated with the 
help of Eq. (|0| ) and Eq. ( p.l2| ). The cosmological parameters were chosen, for illustration, 
as = 10'^, Zdec = 10^, P = —2. The parameter Zi is adjusted in such a manner {zi = 
10^^-^) that the graph goes through the point l{l + l)Ci = 6.4 x 10"^° at / = 10, which 
agrees with observations. Our attention is focused, however, on the oscillations in this 
function. The dashed line shows the same function for the alternative stationary background. 
The cosmological parameters are the same as above, but the calculation is performed with 
the help of Eq. ( |5.8| ) and the coefficients p„, found from the condition ( ^?T| ) of the fair 
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FIG. 2. The solid line depicts the plot of l{l + l)Ci versus I in the physical model, normalized 
such that at I = 10, we have l{l + l)Ci = 6.4 x 10"^'^, which tallies with observations. The (red) 
dashed line is the corresponding plot in the alternative model. Here, we take /? = —2, and the 
redshifts at 772 and rjE to be Zgq = 10000 and Zdec = 1000, respectively. 
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comparison. The remarkable (even if expected) result is that the stationary background of 
gravitational waves does not produce oscillations in the angular power spectrum Ci, whereas 
the non- stationary background does. 

The numerical positions of minima and maxima in the oscillating graph are ordered as 
follows: 

Mimma : h = 137, k = 237, Is = 344, k = 456, k = 569, k = 682, I7 = 796, 



Maxima : li = 161, I2 



269, h 



381, I4 



494, h 



609, L = 723, ly 



839. 



Clearly, these features reflect the oscillations in the metric power spectrum. Judging from 
the mathematical structure of the participating Bessel functions, it is likely that the positions 
of features in the n-space are related to the positions of features in the /-space by a simple 
numerical coefficient of order 1. It is difficult to find out this coefficient analytically, though. 
Remembering that the positions of first features may be displaced by 10-15 percent, as 
compared with the analytical forecast (|3.19| ) , (|3.18| ) , we put the zero-th "would be" maximum 
at 1'^°'^ = 56. Then, our simple analytical formula places the next features in the following 



positions: Z^" = 112, 
448, /^"^ = 504 



168, 



224, / 



max 
2 



280, / 



mm 
3 



336, 



imax 
'3 



392, Z^" 



imin 
1 '5 



560, 



imax 



616, 



672, /^'^^ 



728, = 784, irp 



840. 



Comparing this prediction with the numerically calculated positions, we find them in a fairly 
good agreement. 
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This investigation of gravitational waves provides us with a guidance for the technically 
more complicated case of density perturbations. 



VI. DENSITY PERTURBATIONS AND THE Ci FEATURES 

The general expression for metric perturbations, associated with density perturbations, 
is given by Eq. ( |1.3D with the polarization tensors explained in the Introduction. When 
one is actually writing down the perturbed Einstein equations, it turns out that it is more 
convenient to work with the mode functions h{ri) and /i/(?7) instead of the original mode 

s 

functions hn{v)- The relationship between them is: 

hniv) = \/| (hiri) - hiir])'^ , hn{v) = -^Mv), (6.1) 

where the wave-number index n on the mode functions h{ri) and hi{ri) is implicit. The 

function h{ri) is the purely scalar part of perturbations, it enters hij with the polarization 
1 

structure Pij = 6ij, whereas the function hi{ri) is the purely longitudinal- longitudinal part of 

2 

perturbations, it enters hij with the polarization structure Pij = —niTij/n^. We will follow 
the same strategy as in the case of gravitational waves, and will start from exact solutions 
to the perturbed Einstein equations in different cosmological eras. 



A. Density perturbations at the matter-dominated stage 

The matter-dominated stage is driven by a pressureless matter; possibly, - cold dark 
matter (CDM). The general solution to the perturbed equations at the m stage can be 
simplified by using the available freedom within the class of synchronous coordinate systems. 
By using this freedom, one specializes to the unique coordinate system, which is synchronous 
and comoving with the perturbed pressureless matter. In this coordinate system, the general 
solution is 

h{v) = C^, /^K^) = i^C^i^'(^-^n.)'-^C^2^^^^^, (6.2) 
10 3 (r/ - r]my 

where Ci,C2 are arbitrary complex numbers. The matter density perturbation is 

- = lh,{r^), (6.3) 
and the velocity of matter elements, including perturbations. 



n 

To'-' 



(6.4) 



is by construction zero, i.e.. 
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The above solution is well known since the times of E. M. Lifshitz and can be found in various 
publications, up to possible misprints. For instance, this solution follows from equations 
given in Ref. ^6], if one corrects for a misprint in the last line of Eq. (A8): the second term 



there should actually enter with coefficient 2. 

The adopted choice of the unique coordinate system, which is both comoving and syn- 
chronous, not only greatly simplifies the form of solutions, but is also needed for a proper 
formulation of the 5T/T calculations. As long as the emitter and the receiver are riding 
on the time-like geodesies x* = const, of this perturbed metric, the Sachs- Wolfe integral 
( p.lD is the full answer; there are no extra velocity contributions to this integral. The addi- 
tive velocity contributions arise only if the emitter or the receiver are moving with respect 
to this unique coordinate system, that is, when they are not described by the world-lines 
= const.. 

For the growing solution (namely, the terms with coefficient Ci), the Sachs- Wolfe in- 
tegral can be taken exactly. It appears that astrophysical literature calls by gravitational 
"Sachs- Wolfe effect" only a part of what is actually contained in the Sachs- Wolfe paper . 



Invariably, by the "Sachs- Wolfe effect" are meant only two terms, which are, roughly speak- 
ing, the difference of "gravitational potentials" at the events of emission (E) and reception 
(R). Two other terms in their full formula (43), which are the difference of the "gradients 
of the gravitational potential", are being systematically ignored. Possibly, this happened 
because Sachs and Wolfe addressed one of these gradient terms in the words: "this second 
term is normally small". This second term is indeed small for small wave-numbers, but it 
is in fact dominant for large wave-numbers, which are responsible for the dipole Ci and for 
the Ci multipoles near the peak at / ~ 200. For example, correct implementation of the full 
Sachs- Wolfe formula (43) for calculation of the dipole Ci gives a number that is 5 orders of 
magnitude greater than the number following from the "Sachs- Wolfe effect" counterpart of 
the full formula. The lack of ergodicity on a 2-sphere provides a la uncertainty in the Q's, 



roughly at the level AC/ ~ ^2/ {21 + 1)Q. We say "roughly" because the statistic of the 
underlying random variable is not Gaussian, it is described by the product of an exponent 



and the modified Bessel function Kq ||TJ]. In the case of the dipole Ci, the uncertainty 
amounts to ACi/Ci ~ 0.8. Clearly, this factor-of-two uncertainty cannot cover the 5 orders 
of magnitude disparity in the results; quite simply, the result based on the misinterpreted 



"Sachs- Wolfe effect" is what is wrong. For details, see Ref. . 

The coefficients Ci,C2 in the general solution ( |6.2|) are, so far, arbitrary, but they are 
determined by the previous evolution of density perturbations {d.p.). 

B. Density perturbations at the radiation-dominated stage 

The "master equation" at the e stage is 

zy" + -n^jj = 0, (6.6) 
3 
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where the coefficient 1/3 enters because we have used ci/c= l/-\/3, which is vahd deep in 
the radiation-dominated era. By the time of decouphng, the plasma sound speed decreases 
shghtly below this value, depending on the baryon content, and we will account for this fact 
by returning back to q/c in appropriate places. The general solution to Eq. ( |6.6|) is always 
oscillatory function of time: 



(6.7) 



where -81,-82 are arbitrary complex numbers. All the metric and matter perturbations can 
now be found from solutions (|6.7| ). For metric perturbations, one has: 



and 



h{7]) 



h, = — 
' a' 



rv 

/ udri + Ce 



3h" + 9-h' + n^h 
a 



(6.8) 



(6.9) 



The constant Ce reflects the remaining coordinate freedom at the e stage. The should 
be chosen in such a way that the comoving synchronous coordinate system of the m stage 
joins smoothly to the employed (unique) coordinate system at the e stage; we will discuss 
this specific choice of Ce later on. The constants -81,-82 are still arbitrary and should be 
found from solutions at the i stage. 

The "master equation" at the i stage is 



n — 



0, 



(6.10) 



where 



7 



1 + 



H 

1p' 



and the t-time derivative is related with the ?7-time derivative by cdt = ad?]. For the 
power-law scale factors 0(77) oc l?]]^^^, which we are working with, the function 7 becomes a 
constant, and it drops out of the Eq. (|6.10| ). So, the "master equation" ( |6.10|) is exactly the 
same as equation (2TT) for gravitational waves. By quantum-normalizing the initial metric 
perturbations, and evolving them through the i stage, we finally find that 



-82 = Bd. 



p.- 



(6.11) 



It was shown |TT| that the crucial quantity Bd.p. for density perturbations is related with the 
crucial quantity Bg y,^ for gravitational waves, introduced in Eq. ( |3.9|) , by the relationship 



B, 



d.p. 



Vgb, 



g.w.- 



(6.12) 



In what follows, we will work with B^.p. only and, henceforth, drop the subscript d.p.. 
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Combining all the results together, we write down explicitly the exact solution at the e 
stage, including the required choice of Cg. In doing this, we use the following new notations: 



n 



y2 = -^ir]2-Ve) 



n 



y = 



n 



n 



Zf.q 



(6.13) 
(6.14) 



Then, the exact solution is: 



hi{r]) = 3 A 



Y = -?/2sin?/2 + cos 1/2. 



h{r]) 



A 



[cos y-Y], 



smy fy^cosy y lsin?/2 , 2 

-ay — Y in 1-77 H7rCos?/2 



y Jy y 



Z/2 3 ?/2 3 



(6.15) 
(6.16) 



'1 = -A 



2 2 
— (cos y — F) H — smy — cos y 



y 



y 



(6.17) 



V 



— = -lA — 1= 



-(cosy — Y) + siny 



where 



A = 



inBy/1 + z, 



eq 



(6.18) 



(6.19) 



One can check that all the participating functions, h{ri), h'{ri), hi{rj), h^irf), (5e/eo, f-', join 
continuously with the solution (|6.2| ), (|6.3|) , ( |6.5|) at the transition point rj = rj2. This 
transition fully determines the coefficients Ci and C2: 



"2^ 



sin?/2, C2 



3A r 

5^/2 



(10 - ?,yl) sin 7/2 - IO1/2 cos?/2 



(6.20) 



The oscillatory behavior of Ci, C2, as functions of n, is analogous to the oscillatory behavior 
of the gravitational wave coefficients (|3.12|) and has the same physical origin. The fact that 
Bi —B2 demonstrates that each mode n of the metric perturbations, and the associated 
matter perturbations, at the e stage forms a standing wave pattern. In the limit of short 
waves, ?/ S> 1, one recovers from Eqs. ( |6.17| ), ( |6.18| ) the familiar solutions for standing sound 
waves: 

— Acosy, (6.21) 
Co 



— ^ —lA — ^ 



smy. 



(6.22) 
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C. Perturbations at the last scattering surface 



Having found the quantum-normahzed exact solution at the m stage, we are in a position 
to calculate the metric power spectrum, which is defined by Eq. (|1.71 ). Taking into account 
our mode functions, the spectrum can be written as 



C2 



h--h 
3 



2 



3' " 



(6.23) 



We will calculate this quantity at the last scattering surface 77 = r]E- By that time, the 
second term in the function hiij]) is a factor [(1 + Zdec)/{^ + Zeq)]^^"^ smaller than the first 
term (see Eqs. ( p.2| ) and ( |6.20| )). We neglect this decaying part of the solution, participating 
with the coefficient C2. For the exphcit form of |Cip we use Eq. (|6.2CI|) and Eq. ( [6.191) . 
Then, we obtain 

2. . n^\B\\l + z,,) ^sini/2 V (300 - 20p^t/| + A|) 
^ ^"'^-^ = 2^ 48l^ [-^) 200 ' ^'-''^ 

where we have introduced the quantity 



P 



2x/3^rT 



eq 



related to a similar quantity, 6, from the g.w. case, hy p = \^b. The spectrum certainly 
retains its primordial form in the band of long waves n <^ a/1 + Zdec- Taking into account 
Eq. ( |6.12| ) and the numerical value of C = y/2A7ilpi for density perturbations, one finds 
that the primordial spectrum of metric perturbations associated with density perturbations, 
Eq. ( |6.24| ), is a factor 9/16 lower than its gravitational wave counterpart, Eq. ( ^.15| ). In 
particular, for (3 = —2, one finds 

h\n, r]E) ^ n < ridec = Vl + Zdec- (6.25) 

IT Ijj [1 + Zeq) 

For relatively short waves, n/n^ ^ 1, the crucial part of the power spectrum ( p.24| ) is 
the modulating (transfer) function 

n\ /sint/2\^ sin^ (n/rZc) 



The primordial metric spectrum is encoded in the factor ?T,^|i?p. Whatever this spectrum 
is, the modulating function leaves it intact at large scales, but bends the spectrum down 
and introduces oscillations at smaller scales. In Fig. ^, we show the metric power spectrum 
h'^{n,riE), up to numerical coefficients. (To avoid confusion, we emphasize again that this is 
the spectrum of the (squeezed) metric perturbations associated with density perturbations, 



and not the gravitational-wave spectrum ||23[.) Specifically, by a solid line, we plot the 
function p{x,p), where 
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smx 



X 



300 - 2Qp^x^ + p^x 



,4 4 



X = Di = njrLc and, for illustration, we take p = 8. By a dashed line, we plot the function 
M^(x), multiplied by the artificial numerical factor 10® in order to facilitate the visual 
comparison of maxima and zeros in the two graphs. 



FIG. 3. The plot depicted by the solid line is that of 
shows the behavior of M'^(x) x 10®. 



versus x. The dashed line 




We now turn to the ordinary matter perturbations at the last scattering surface. The 
photon-electron-baryon fluid is gravitationally subdominant at t] = rjE- The fluid does not 
significantly contribute to metric perturbations, but it retains its own perturbations. The 
plasma speed of sound is given by 

^ (6.27) 



c ^3(1 + R) 

where R = 3pb/4:p^ ^ 27VLbh? [ITS, [EDI. For the popular value Vtf.h'^ ^ 0.02, it means that 



ci/c decreases from the nominal value 0.58 to approximately 0.47. The plasma standing 
waves ( |6.21D , (|6.22|) , continued to the decoupling era r] = tie, take the form: 



^^"f^Acos-, (6.28) 



— ^ -iA-— sin — , (6.29) 

c c n n. 
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where 



ris = — VI + Zdec = —Udec- (6.30) 

Cl Cl 

[The velocity is always defined with respect to the unique synchronous coordinate system, 
which is comoving with the gravitationally dominant pressureless (dark) matter.] At the 
same time, the leading metric perturbation is given by 

hi{vE) = -^A^±^-sm-. (6.31) 

5 1 + Zdec rie Tic 

There are a number of differences between the metric perturbations and the plasma 
perturbations at the last scattering surface. First, the amplitude of hi{r]E) is, at least, 
a factor 3(1 + Zeq)/[5{1 + Zdec)] greater than the amplitudes of (5e^/e^ and /c, near the 
most interesting scales n ^ Uc- After all, the original motivation for the introduction of a 
cosmological dark matter was precisely this: to avoid conflicts with 6T/T observations by 
allowing the plasma perturbations at decoupling to be small, but, nevertheless, to be able 
to develop the large scale structure of luminous matter, at the expense of large gravitational 
fleld perturbations driven by the dark matter. So, we have to pay the price for this idea 
by exploring in more detail the consequences of large metric perturbations for the CMB 
anisotropies. 

Second, the characteristic frequencies Uc and Us are different. Their ratio is 



Uc 




(6.32) 



The Zeq is given by 1 + ~ 4 x 10 fi^/i |TP|. For the popular values Vt^ = 0.3, h = 0.7, 
this amounts to ~ 6 x 10^. So, the ratio rtc/ns can be a number close to 4. 

Third, although the sound waves before the decoupling are standing waves, they are 
still not processed by the quick drop of the sound speed to zero at the decoupling. This 
processing will later lead to the baryonic matter power spectrum modulations known as 
the Sakharov oscillations They would have taken place even in laboratory conditions. 



where gravity plays no role. The Sakharov oscillations are important for the formation of 
oscillating features in the luminous matter power spectrum, but they are unlikely to be 
directly responsible for the peaks and dips in the observed C;'s. In a broad sense, the 
periodicity in the metric power spectrum, related to the transition t] = ri2, can also be called 
Sakharov oscillations, but this is not what was originally meant by the Sakharov oscillations. 
In short, the zeros in the metric power spectrum are "frozen" zeros, they are determined by 
M'^in/ric)'-, whereas the zeros in the plasma power spectrum, at the times before decoupling, 
are still "moving" zeros; they change their positions at slightly different moments of time 
T] = const. [| 



The notion of the "moving zeros" was suggested by J. Peebles in a private correspondence, May 
1990. 
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Fourth, the wave-number periodicity in the metric power spectrum is governed by the 
sine function, whereas the periodicity in Eq. ( |6.28 ) is governed by the cosine function. 
Presently, there exists a tendency to distinguish between the "acoustic peaks" in the Ci 
(supposedly caused by Eq. (|6.28|) and by the "effective temperature") and the "Doppler 
peaks" (supposedly caused by the velocity in Eq. ( |6.29| )). The authors of |^ emphasize that 
"the acoustic peaks are not "Doppler peaks" " , arguing that the irrotational velocity cannot 
produce strong peak structures in the Ci spectrum. They say that "the observed peaks must 
be acoustic peaks" and they give the ratio of the peak locations: £i : £2 : ^3 ~ 1 : 2 : 3 . 
So, the main contenders for the explanation of the peaks seem to be the sine function in the 
metric power spectrum and the cosine function in the "acoustic peaks" . 

Before proceeding to the discussion of peaks and dips, we need to make one more 
comment. It was shown above that the primordial power spectra of gravitational waves 
and density perturbations are of the same order of magnitude, with some small numer- 
ical preference for gravitational waves. In particular, this is true for the flat spectra 
{j3 = —2), as demonstrated in Eq. (p.l6|) and Eq. (|6.25|) . Therefore, the lower order 
CMBR anisotropies (starting from the quadrupole moment C2) are expected to be of the 
same order of magnitude ||Tl|. One should be aware that the story is dramatically differ- 
ent in inflationary scenario. The "standard result" of inflationary scenario [PS|-|5B[], pre- 
dicts the infinitely large density perturbations, in the limit of the flat spectrum (that is, 
the Harrison-Zel'dovich-Peebles spectrum, with spectral index n = 1, parameter jS = —2, 
and the relationship between them being n = 2(3 + 5), through the set of evaluations: 

composing the ratio of the gravi- 
tational wave amplitude Ht to the predicted divergent amplitude of the scalar metric per- 
turbations hs (the so called "consistency relation": hr/hs ~ \/l — n), inflationary theorists 
substitute their prediction of arbitrarily large density perturbations for the claim that it is 
the amount of gravitational waves that should be zero, or almost zero, at cosmological scales 
and, hence, down to laboratory scales. This claim has led to many years of mistreatment 
of a possible g.w. contribution to the CMBR data. It is only in a few recent papers (for 
example, ^7j) that the inflationary "consistency relation" is not being used when analyzing 



the CMBR and large scale structure observations, with some interesting conclusions. For 
the latest statement that the initial spectrum of gravitational waves is "constrained to be 
small compared with the initial density spectrum" see the latest article praising inflationary 
predictions (for instance, pSf). For the critical analysis of the "standard inflationary result" 
see the end of Sec. |V1| in and references therein. 



D. Peaks and dips in the angular power spectrum 

We will now analyze the zeros and maxima of the metric power spectrum P{x,p) shown 
in Fig. 1^. The crucial periodic dependence is provided by the function sin^(a;) = sin^(n/nc). 
We will use this function for our analytical evaluation, in full analogy with the case of 
gravitational waves. The positions of maxima and zeros are determined by the rules: 

Maxima : x^'' = f^^^ + = a;^"^(2fc + 1), A; = 0, 1, 2, 3, (6.33) 
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Zeros : x^" = ^2fc, x™'' = xf^'k, k = 1,2, 3, .... (6.34) 

Obviously, the zeros of the function sin^(x) are exactly the same as the zeros of the full 
function p{x,p). But the positions of maxima are somewhat different. The difference is 
significant for the zero-th maximum, but it fully disappears for later maxima. The locations 



of the first few maxima, derived from the simple analytical formula (|6.33|) , are: x™""^ 



1.57, x™"^ = 4.71, x^°^ = 7.85, x^"^ = 11.00. At the same time, accurate positions from 
the numerical calculation are: x[}^"^ = 2.05, x^'^^ = 4.92, x^"^ = 7.98, x^'^^ = 11.09. Thus, 
formula ( |6.33| ) predicts the positions of the first two maxima somewhat to the left (smaller 
x) than they should actually appear, but the positions of zeros and further maxima are 
described very well. In terms of the percentage corrections, the zero-th maximum, derived 
from ( |6.33|) , should be shifted to the right by 30 percent, and the first maximum should be 
shifted to the right by 4 percent. 

Accepting Zeq = & ^ 10^, one obtains n^. = 268. With this Uc and x[J^°^ = 1.57, the 
position of the zero-th maximum in the n space would be, according to ( 6.33| ), at n"^"-^ = 421. 



Positions of all the subsequent features in the power spectrum follow from the general rules 
( |6.34| ) and (|6.33|) . The problem now is to relate these features in the metric power spectrum 
with the peaks and dips in the angular power spectrum l{l + l)Ci. Judging grom the 
previous numerical experience [0, the characteristic features of the metric power spectrum 



are reflected in the /-space via a numerical coefficient a close to 1/2: / = an. Accepting 
the provisional value a ~ 1/2, the location of the zero-th peak in the / space would be near 
/o = 210. This is a satisfactory intermediate result, but we want to do better. Remembering 
that the position of the zero-th peak, following from the analytical formula ( |6.33|) , should 



be shifted to the larger values of /, we place our zero-th peak at = 170. The 30 percent 
correction of this number shifts the zero-th peak to Z™"^^ = 221. Of course, we keep an eye 
on the actually detected peak in this region. Our aim is to derive the full structure of peaks 
and dips in the angular power spectrum l{l + l)Ci from the simple analytical formulas ( |6.33 ) 



and (|6.34| ), allowing only for the 30 percent correction to the zero-th peak and the 4 percent 
correction to the first peak. Following this strategy, we formulate our full forecast: 

Peaks : = 170(221), /r"" = 510(530), Z^^' = 850, Zg"^"^' = 1190, /r"" = 1530 , (6.35a) 
Dips : Z™" = 340, /r" = 680, Z^" = 1020, /f " = 1360, /™" = 1700 . (6.35b) 



As a consequence of Eqs. (|6.33|) and (|6.34| ), the general rule for the peak positions is 



1 : 3 : 5 : 7..., for the dip positions: 1:2:3:4..., and the dips appear between the peaks at 

imin 1 f jmax i ]max\ 

Everywhere in this paper, both for gravitational waves and density perturbations, we 
perform calculations under the simplifying assumption that the Universe is spatially flat. 
It is obvious, however, that neither the generating mechanism itself nor the results, for 
wavelengths comfortably shorter than the putative curvature radius, depend on this simpli- 
fication. The unaccounted factors, such as the possible presence of a spatial curvature, or 
a A term, or a "quintessence" , or a "dark energy" , can move the entire structure of peaks 
and dips, but these factors can hardly change the general rules for their relative positions. 
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One should note that what is following from our classification as the "zero-th gravita- 
tional peak", which we place at Z™""^ = 170 plus the correction shifting it to /™"^ = 221, 
is often interpreted as the "first Doppler peak" or the "first acoustic peak" . The notion of 
the "zeroth Doppler peak" was introduced and discussed by Weinberg [|l9l, pOf. In general. 



all three sources - gravitational field perturbations, intrinsic temperature variations, and 
velocities - contribute to the peak structure. The gravitational field contribution is repre- 
sented by the Sachs- Wolfe integral ( p.l[ ), while the two other sources are represented by Eqs. 
( |6.28| ), ( |6.29|) . However, the raising function Z(Z + 1)C/ would not have turned down with- 
out the modulating function M'^in/ric) [0, so we focus our attention on the gravitational 
contribution. 

The numerical graph of Fig. |^ shows also a little depression at Xdep = 0.41. This depres- 
sion arises entirely due to the polynomial term in p{x,p) rather than from the modulating 
function M^(x). Accepting the same value ric = 268, this feature corresponds to n^ep = 110. 
Assigning some significance to this feature, and following the same logic as before, we have 
to conclude that this depression in the metric power spectrum may be reflected as a small 
local minimum in the angular power spectrum. Applying the numerical factor a = 1/2, this 
minimum is expected to be seen around Idep ~ 55. This may be one of the areas in the I 
space to analyse closely in the future experiments, such as MAP and PLANCK. 

To compare our forecast with observations, we take for the face value the central positions 
of peaks and dips reported by de Bernardis et al. ||l^. We take the liberty of calling their 
Peak 1 as our zero-th peak. Peak 2 as the first peak, and so on. The reported measured 
positions follows: 

Peaks : Z™"^' = 213, /r"" = 541, /r"" = 845 , (6.36a) 
Dips: /™" = 416, /™" = 750. (6.36b) 

Their forecast for the next features is as follows: 

Peaks : /^"^ = 1139, /f''^ = 1442 , (6.37a) 
Dips : /^^'^ = 1025, " = 1328, = 1661 . (6.37b) 

Comparing the observed positions ( |6.36| ) with our formulas (|6.35| ), we find them in reasonably 
good agreement. The peaks and dips appear, at least roughly, in the right positions. On the 
other hand, the periodic function cos^{n/ns), appropriate for the "acoustic peaks", implies 
the reversed rules for the dip and peak positions: 

Dips: /r" = C"(2A; + 1), A; = 0,1,2,3,..; Peaks: Zr"" = ^r'''^, A; = 1, 2, 3, .... 

So, the structure is supposed to start from the zero-th dip, the dip positions are ordered 
as 1 : 3 : 5 : 7..., the peak positions are ordered as 1 : 2 : 3 : 4..., and the peaks appear 
between neighboring dips at Z™"^ = (l/2)(/p" + /^'f). If the first acoustic peak is at / 213, 
the second one is supposed to be at / ~ 426, almost in the same place where observations 
indicate the first dip. Most importantly, there is no observational evidence whatsoever for 
the zero-th dip. We do not see how the observed structure ( |6.36| ), can be explained by the 
acoustic peaks. 
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The forecast ( |6.35| ) on one side, and the forecast (|6.37|) on the other side, go out of phase 
at late features. We place our fourth peak inbetween the positions where de Bernardis et al. 
I7| place their and If these features are not washed out by damping p!9| , |20|j28 



the MAP and PLANCK missions will provide the answer. So far, we tentatively conclude 
that the structures in the angular power spectrum are caused by squeezing in the primordial 
gravitational field perturbations associated with the density perturbations. 
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Note Added in Proofs: 

The latest CBI observations (T.J.Pearson et al., ^stro-ph/0205388] ) have detected four 
peaks, at / ~ 550, 800, 1150, 1500, and four dips, at / ~ 400, 700, 1050, 1400. These positions 
are in a very good agreement with the theoretical formula (6.35) of the present paper. We 
interpret this data as confirmation of our conclusion that it is gravity, and not acoustics, 
that is responsible for the observed structure. 
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